Abnormal effective connectivity in migraine with aura under photic stimulation 
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Migraine patients with aura show a pecuhar pattern of visual reactivity compared with those of 
migraine patients without aura: an increased effective connectivity, connected to a reduced syn- 
chronization among EEG channels, for frequencies in the beta band. The effective connectivity is 
evaluated in terms of the Granger causality. This anomalous response to visual stimuli may play a 
crucial role in the progression of spreading depression and clinical evidences of aura symptoms. 

PACS numbers: 87.19.L-, 05.45.Tp, 05.45.Xt, 42.66. Lc 
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Important information on the structure of complex 
systems can be obtained by measuring to what extent 
the individual components exchange information among 
each other. Transfer entropy and Granger causality 
have emerged in recent years as leading statistical 
techniques to detect cause-effect relationships between 
time series; they are equivalent in the case of Gaussian 
stochastic variables Q. These approaches provide fur- 
ther insights on the architecture of complex systems, in 
addition to those from correlation and synchronization 
analysis (8l-[l0|. In neuroscience, interdependencies esti- 
mated by correlation or spectral coherence are referred 
to as functional connectivity, whilst effective connectiv- 
ity is a notion related to Granger causality or transfer 
entropy [Tlj : the relationship between functional and ef- 
fective connectivity in the cortex represents a significant 
challenge to present-day neuroscience [T2| - [l5| . 

Migraine, an incapacitating disorder of neurovascu- 
lar origin, consists of attacks of headache, accompanied 
by autonomic and possibly neurological symptoms [l6| . 
This pathology affects a relevant fraction of the general 
population and represents a social problem. The study 
of phase synchronization in EEG rhythms showed a pat- 
tern of alpha rhythm (8-12.5 Hz) hyper-synchronization 
under repetitive flash stimulation in migraine without 
aura patients, opposite to a de-synchronization trend in 
non-migraine subjects 117]; this nonlinear EEG pattern 
was found to be modulated by anti-epileptic drugs [l8| . 
Approximately one-third of people who suffer migraine 
headaches perceive an aura (visual abnormality lasting 
10-30 minutes) as a sign that the migraine will soon oc- 
cur [l^. There are evidences [13] of relationships be- 
tween migraine aura and the spreading depression (SD) 
phenomenon (a wave of electrophysiological hyperactiv- 
ity followed by a wave of inhibition ^21*] ) but the link be- 
tween SD and headache is far to be explained in humans 
[2^ . SD in clinically less conspicuous or extended brain 
regions may be the trigger of mig raine attacks ostensibly 
without a "perceived" aura [23|, so it may be plausi- 
ble that, migraine with and without aura may differ for 
those neuronal factors favoring the clinical expression of 
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FIG. 1: (Top left) The hnear Granger causahty of EEGs 
filtered in the beta band, averaged over pairs of channels 
and over subjects in the classes, is depicted for aura patients 
(stars), no aura patients (empty circles) and controls (empty 
squares) for basal condition and as a function of the frequency 
of stimulations. We also make the supervised analysis (hy- 
pothesis testing) of the linear Granger causality, averaged 
over pairs of channels, in the beta band. The probabilities 
that the measured values from two classes were drawn from 
the same distribution, evaluated by t test, are depicted: aura 
vs controls (top right) aura vs no aura (bottom left) and aura 
vs controls (bottom right). Dashed lines correspond to the 
significance threshold at 5% after Bonferroni's correction. 



SD phenomenon. 

The question we address here is the following: how 
does the response of aura migraine patients, to external 
stimuli, differ from those of patients without aura? We 
show that, in presence of visual stimulations, migraineurs 
with aura show a pattern of increased effective connec- 
tivity between EEG channels in the beta band, and cor- 
respondingly a decrease of the synchronization between 
channels. This variation of connectivity, due to stimuli, 
is statistically significant for migraine with aura: we re- 
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FIG. 2: (Top left) The nonlinear Granger causality of EEGs 
filtered in the beta band, averaged over pairs of channels 
and over subjects in the classes, is depicted for aura patients 
(stars), no aura patients (empty circles) and controls (empty 
squares) for basal condition and as a function of the frequency 
of stimulations. A Gaussian kernel with a = 10 is used. We 
also make the supervised analysis (hypothesis testing) of the 
nonlinear Granger causality, averaged over pairs of channels, 
in the beta band. The probabilities that the measured val- 
ues from two classes were drawn from the same distribution, 
evaluated by t test, are depicted: aura vs controls (top right) 
aura vs no aura (bottom left) and aura vs controls (bottom 
right). Dashed lines correspond to the significance threshold 
at 5% after Bonferroni's correction. The ANOVA tak- 
ing the nonlinear Granger causality as variable, the diagnosis 
(migraine with aura vs migraine without aura vs controls) 
and the frequency of stimulation (9-18-21-24-27 Hz) as fac- 
tors, yields the F values 377.35, 8.36 and 3.38 for diagnosis, 
frequency and interaction respectively; all the three factors 
are thus recognized as significant. 



mark that in our knowledge neurophysiological patterns 
separating migraine with and without aura have been 
rarely detected 12411 . nor relevant differences in regard to 
visual reactivity |25| . 

Our data are as follows. EEG is recorded from 19 
patients (7 males, 20-44 age) affected by migraine with 
aura, 19 (4 males, 21-45 age) patients affected by mi- 
graine without aura, and from 11 healthy subjects (con- 
trol group, 3 males, 20-46 age). All patients are in the 
interictal state, the time from the end of the last attack 
being at least 72 h. No patient was under preventive 
treatment nor had assumed symptomatic drugs in the 
72 hours preceding the recording session. During the 
acquisition, flash stimuli are presented at a rate of 9-18- 
21-24-27 Hz; also EEG in the absence of stimuh (base) 
is recorded. Each frequency of stimulation is delivered 
by a flash with 0.2 J luminance for about 20 sec. EEG 
data are recorded by six scalp electrodes: two occipital 
channels (01 and 02), two parietal ones (P3 and P4), a 
central electrode (Cz) and a frontal one (Fz); the sam- 
pling rate is 256 Hz, and the EEG is digitally filtered off 



FIG. 3: For each electrode, the incoming (top) and outgoing 
(bottom) Gaussian kernel Granger causality in the beta band, 
corresponding to 18 Hz stimuli and averaged over subjects in 
classes, is depicted; aura patients (stars), no aura patients 
(empty circles) and controls (empty squares). 

line by a filter with a band-pass 0.3-30 Hz. 

Next, we describe our findings. Firstly we investi gate 
the alpha band hyper-synchronization phenomenon [l7[ 
in presence of flash stimuli. The pattern of [l3| is con- 
firmed for migraineurs without aura, whilst patients with 
aura do not show alpha band hyper-synchronization in 
presence of light stimuli: in this range of frequencies, 
aura patients seem to behave as controls. 

Moving to higher frequencies (beta band, 12.5-30 Hz), 
instead, we find a peculiar pattern of visual reactivity for 
aura patients. We set the order of the regression model 
equal to TO = 6 ^26'| and evaluate the linear Granger 
causality in the beta band [l^l : linear causalities are weak 
and show significant differences among classes only at 18 
Hz stimulations, see figure ((IJ. On the other hand, evalu- 
ating the nonlinear Granger causality among the filtered 
EEG signals with the kernel approach described in Q 
and a Gaussian kernel, we find that aura migraine pa- 
tients exhibit increased values of causality in presence of 
stimuli, whilst controls and no aura patients do not show 
significant variation w.r.t. basal conditions. In particular 
the discrimination between aura and no-aura patients is 
excellent in presence of flash stimuli. Figure refers to 
cr = 10 as the width of the Gaussian kernel, however the 
results are robust to variations of a. 

A topographic analysis is also performed: for each elec- 
trode, we evaluate the total incoming causality (the sum 
of the causalities from the other electrodes to the elec- 
trode under consideration) as well as the total outgoing 
causality (the sum of the causalities from the electrode 
under consideration to the other electrodes). In figure 
([3]) we depict these quantities for the three classes at 18 
Hz stimulations in the beta band: it shows that the phe- 
nomenon is diffuse over the scalp. 

Use of the Gaussian kernel ensures that all orders of 
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FIG. 4: (Top left) The nonlinear Granger causality of EEGs 
filtered in the beta band, averaged over pairs of channels 
and over subjects in the classes, is depicted for aura patients 
(stars), no aura patients (empty circles) and controls (empty 
squares) for basal condition and as a function of the frequency 
of stimulations. The polynomial kernel with p = 2 is used. We 
also make the supervised analysis (hypothesis testing) of the 
polynomial Granger causality, averaged over pairs of channels, 
in the beta band. The probabilities that the measured val- 
ues from two classes were drawn from the same distribution, 
evaluated by t test, are depicted: aura vs controls (top right) 
aura vs no aura (bottom left) and aura vs controls (bottom 
right). Dashed lines correspond to the significance threshold 
at 5% after Bonferroni's correction. 



nonlinearities are taken into account. If one considers 
only nonlinearities up to the second order, weaker causal- 
ities are detected and less discriminating power is ob- 
tained, see figure ^ where we describe the application 
of nonlinear Granger causality with a polynomial ker- 
nel of degree two. We conclude that a relevant amount 
of nonlinear information transmission characterizes this 
phenomenon. 

Turning to synchronization, we consider the Pearson 
linear correlation between channels and find that the in- 
creased flow of information, due to flash stimuli, is con- 
nected to weakening of the correlations between them. 
In flgure ([5]) we depict the linear correlation between sig- 
nals filtered in the beta band, averaged over all pairs of 
channels; in presence of light stimulations the strength of 
correlations decreases for aura patients. We find similar 
results also in terms of the coherence function averaged 
in the beta band, see figure dBJ, as well as for the beta 
band phase synchronization 29]. 

We also quantify the separation among classes in terms 
of the ROC area [30| , which is directly related to the sepa- 
ration of two conditional distributions, and measures the 
discrimination ability of the forecast. For any frequency 



of stimulation, using the Gaussian kernel Granger causal- 
ity we obtain a roc area equal to 0.87 for aura vs no-aura 
patients. A similar result is found using the linear cor- 
relation, roc area equal to 0.82 for aura vs no-aura at all 
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FIG. 5: (Top left) The linear correlation of EEGs fiftered in 
the beta band, averaged over pairs of channels and over sub- 
jects in the classes, is depicted for aura patients (stars), no 
aura patients (empty circles) and controls (empty squares) for 
basal condition and as a function of the frequency of stimula- 
tions. We also make the supervised analysis (hypothesis test- 
ing) of the linear correlations, averaged over pairs of channels, 
in the beta band. The probabilities that the measured val- 
ues from two classes were drawn from the same distribution, 
evaluated by t test, are depicted: aura vs controls (top right) 
aura vs no aura (bottom left) and aura vs controls (bottom 
right). Dashed lines correspond to the significance threshold 
at 5% after Bonferroni's correction. The ANOVA [11], tak- 
ing the linear correlation as variable, the diagnosis (migraine 
with aura vs migraine without aura vs controls) and the fre- 
quency of stimulation (9-18-21-24-27 Hz) as factors, yields the 
F values 143.42, 3.67 and 0.74 for diagnosis, frequency and in- 
teraction respectively; diagnosis and frequency of stimulation 
are recognized as significant factors. 



frequencies. 

Summarizing, we have described for the first time a 
neurophysiological pattern which seems peculiar of mi- 
graine patients perceiving visual aura, where the loss of 
synchronization between channels, due to light, induces 
stronger statistical causal connections among them, dif- 
fuse over the scalp. This pattern is characterized by non- 
linear transfer of information among channels and pro- 
vides excellent discrimination between aura and no-aura 
patients. The biological implications of this complex phe- 
nomenon in facilitating SD progression and aura symp- 
toms perception is the challenge for a better understand- 
ing of migraine pathophysiology. 
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